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Compact astrophysical objects are a window for the study of strongly interacting nuclear matter 
given the conditions in their interiors, which are not reproduced in a laboratory environment. Much 
has been debated about their composition with possibilities ranging from a simple mixture of mostly 
protons and neutrons to deconfined quark matter. Recent observations on the mass of two pulsars, 
PSR J1614-2230 and PSR J0348+0432, have posed a great restriction on their composition, since 
the equation of state must be hard enough to support masses of about at least two solar masses. 
The onset of quarks tends to soften the equation of state, but it can get substantially stiffer since in 
the high-dense medium a repulsive vector interaction channel is opened. Nevertheless, we show that 
once gluon effects are considered, the equation of state of strange stars formed by quark matter in the 
color-flavor-locked (CFL) phase of color superconductivity becomes softer decreasing the maximum 
stellar mass that can be reached. This may indicate that strange stars made entirely of CFL matter 
can only be favored if other interactions, as the one corresponding to the vector channel, are taken 
into consideration and are large enough. 

PACS numbers: 21.65.Qr, 26.60.Kp, 97.60.Jd 


I. INTRODUCTION 


Recent very precise mass measurements for two compact objects, PSR J1614-2230 and PSR J0348+0432 with 
M = 1.97 ± 0.04Mq [lj and M = 2.01 ± O.O4M 0 j2|, respectively, where M 0 is the solar mass, have provided an 
important and reliable way to constrain the interior composition of neutron stars. These high mass values imply that 
the equation of state (EOS) of the corresponding stellar medium should be rather stiff at high densities. 

On the other hand, in the highly dense cores of compact objects, the neutron-rich matter can give rise to more 
degrees of freedom, like hyperons, and perhaps even transitioning to a quark matter phase (see [3] for review). Given 
that cold strange quark matter has been argued to be absolutely stable , a phase transition could occur in the core 
of a compact star that would quickly favor quark matter in all its interior, thus giving rise to a strange star. 

Although it has been known for more than a decade that perturbative QCD predicts maximum stellar masses for 
strange stars larger than 2M 0 II , there were claims [3, H|, already from the first indications of the existence of 
2 Mq stars Q, that quark matter might have to be ruled out as a possible core phase. One reason was that the density 
of compact stars is not expected to be sufficiently high to validate a perturbative approach of the strong interaction. 
On the other hand, nonperturbative calculations based on simple QCD phenomenological models like the NJL with a 
four-fermion quark-antiquark channel, were found to render a too soft EOS, incapable to stabilize a 2M 0 star against 
the gravitational collapse (see for example 0)- However, when other interactions, which are always present in a 
dense medium of quarks [ll|, such as those corresponding to the diquark channel [l2| and/or the vector channel, 
are taken into account, the EOS can become stiff enough. Because of these findings, quark matter was back in the 
competition as a possible core phase of massive stars (see for instance El)- One may wonder if this last conclusion 
could be challenged in turn by other effects that might be part of the complicated physics of super dense quark matter. 

Both, weakly coupled QCD and NJL models inspired on the one-gluon exchange interaction of QCD have predicted 
that the most favored phase of color superconductivity at asymptotically large densities is the color-flavor-locked (CFL) 
phase 0 - In the case of the NJL approach, the lack of confinement of this model is reflected in the introduction 
of an undetermined vacuum contribution to the EOS that enters as an independent parameter in the form of a 
bag constant. Moreover, in this approach, gluons degrees of freedom are usually disregarded as negligible at zero 
temperature. Nevertheless, as discussed in the present paper, gluons can affect the EOS of the CFL phase. The 
reason is connected to the fact that the polarization effects of quarks in the CFL background embed all the gluons 
with Debye (mu) and Meissner (tom) masses that depend on the chemical potential /i [15l |. 
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As a consequence, the gluons in the CFL phase acquire nonzero rest energies that yield to a positive contribution to 
the system’s energy density and thus to a negative contribution to the pressure, thereby softening the EOS. In (JT]) , 
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Nf is the number of flavors of massless quarks, and g the quark-gluon gauge coupling constant. 

The main goal of this paper is to investigate the gluon’s effects on the EOS and mass-radius (M-R) relationship of 
quark stars with CFL color superconductivity using a gauged NJL model. We will show that the contribution of the 
gluon rest energy to the CFL thermodynamic potential softens the EOS and requires a relatively high value of the 
vector interaction if the observed 2M 0 masses for this class of compact objects is to be reached. 

The paper is organized as follows: In Sec. II, the gauged-NJL model is introduced and its one-loop thermodynamic 
potential in the liard-loop approximation is found. In Sec. Ill, we investigate the EOS of quark matter in the CFL 
phase, with and without gluon contributions, for different values of the vector interaction coupling. Then, in Sec. IV, 
the M-R relationships for compact stars having quarks in the CFL phase, with and without gluon contributions, are 
graphically obtained for a suitable range of parameters and the implications of the results are discussed. In Sec. V, 
the main conclusions derived from the paper outcomes are listed. 


II. GAUGED-NJL MODEL WITH VECTOR INTERACTIONS 


Consider the three-flavor gauged-NJL Lagrangian at finite baryon density, which is invariant under SU C ( 3) x 
SUl( 3) x SUu(3) x t/y(1) symmetry 
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In ([2]). the quark fields ipf have flavor (1 = u,d,s) and color (a = r,g,b ) indexes. We consider the quark-antiquark 
channel with coupling constant Gs, the ’t Hooft determinant term that excludes the Ua{ 1) symmetry jT(| with 
coupling constant K , the diquark channel P ^ = C^e^e^rj with coupling Gd, and the repulsive vector channel 
with coupling Gy. This last interaction naturally appears after a Fierz transformation of a point-like four-fermion 
interaction with the Lorentz symmetry broken by the finite density 0- We neglect the current quark masses mo?; since 
our main interest is to explore the stable CFL phase, which can only exist in the density region where m 2 < 2pA cfl- 
In this way, we avoid the chromomagnetic instabilities of the gCFL phase [l8|. Thus, our results are valid in the 
density domain where the contribution p 2 Aq FL to the EOS is more important than that of the s-quark mass. Hence, 
it is consistent to neglect m s in our calculations. To open the door to possible m s effects relevant at lower densities, 
the stable ground state that removes the chromomagnetic instability at those densities would have to be determined 
first, a problem out of the scope of this paper. Despite several propositions |l9j]-[2l|, there is still not a definitive 
accepted answer. 

The gluon Lagrangian density £ G is 
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with G A v = d g G A - d u G A + gf ABC G B G ^, the gluon strength tensor; £ gauge, a gauge-fixing term; and £ gho st = 
—r] A ^d ,i {d tl ri A + gf ABC G^r] C ), the ghost-field Lagrangian. The coupling between gluons and quarks occurs through 
the covariant derivative 


D^=d^- igT A G A , (4) 

where T A = X A /2 7 A = 1 — 8, are the generators of the SU C { 3) group in the fundamental representation (X A are 
the Gell-Mann matrices). The model © is nonrenormalizable due to the four-fermion interaction terms, so a cutoff 
A needs to be introduced in all the calculations to regularize the theory in the ultraviolet region. The parameter A 
defines the energy scale below which this effective theory is valid. 

We are interested in exploring the region of densities large enough to have a CFL phase. Then, we can use that 
the chiral condensate has already vanished and hence consider only the expectation values for the diquark condensate 
A, ; = {ip T Prjif)) and the baryon charge density p = {'ijj'yoip). One can now bozonize the four-fermion interaction 
via a Hubbard-Stratonovich transformation and then take the mean-field approximation to obtain the mean-field 
Lagrangian 
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where fi = /z — 2 Gyp. 

The partition function in the mean-field approximation is then given by 

Z M f = N J V4>Vil)VGVrfVrtef CMFdix (6) 

Even though the mean-field Lagrangian fl5j is quadratic in the fermion fields, the presence of fermion-fermion terms 
does not allow to use the conventional formula to integrate in the fermion fields. However, we can use the well-known 
Nambu-Gorkov (NG) method to overcome this difficulty. For that, we rewrite Cmf in terms of ip and its charge 
conjugate ipc = Cip T , with C = * 7 2 7 o the charge operator, and introduce NG spinors 

Hi 

Now we can transform to Euclidean variables and consider a space with volume V and Euclidean time of length 
/? = 1/T, to obtain the finite-temperature, mean-field partition function 
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the inverse NG fermion propagator in the color superconducting medium. Notice that Sq 1 [x) depends on the gluon 
field through the covariant derivative. The symbol = A V P V denotes the CFL gap matrix and = 7o[$ + ]^7o 
[22l ]. Replacing in © the covariant derivative by the derivative implies Sq 1 (a;) —> 5 ( 7 1 (x)|g = o, where 5 q 1 (o;)|g=o is 
the mean-field inverse fermion propagator of the CFL phase of the non-gauged NJL model. 

Performing the integral in the NG spinors in © we obtain 


Z M f = e~ 0vnMF = Afe~ 0V( ~~^ !L+Gvp2) J VG*V'l ] Vr 1 e- T{G)+ S c ° dix , (10) 

where 
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The effective action T(G) can be expanded in powers of the gluon fields as 

OO 
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where the tensors n^A 2 ’;;^"- are n-vertex (with vertex gT^y^) one-loop polarization operators with internal lines 
of fermions with propagator 
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Here 

5 ± = ([5 0 ± ]- 1 -S ± )- 1 , (14) 

denotes the propagators for quasiparticles and charge-conjugate quasiparticles, with Sq the propagator for massless 
free quarks, and 
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the quark self-energy generated by the exchange of particles or conjugate particles with the color superconducting 
condensate. The off-diagonal elements of (1131) are given by 

S ± = -S'^$ ± S ,± . (16) 

Except for the replacement of p by p. due to the vector interaction, the zero-order term Slo in (112p exactly reproduces 
the fermion determinant part of the conventional (no gluons and no vector interaction) mean-field thermodynamic 
potential of the CFL-phase [23| . Going to momentum space and summing in the Matsubara frequencies, the quark 
contribution to the thermodynamic potential can be written as 

Q A 2 1 /-A _ o 1A 2 

^cfl{T) = D 0 + — Gvp 2 = J dpp 2 yie,| + — ln(l + e — Gyp 2 , (17) 

where the e^’s are the dispersion relations of the quasi-quarks/antiquarks in the CFL phase and we used that in the 
chiral limit Ai = A 2 = A 3 = A. 

In the zero-temperature limit the logarithm vanishes and we obtain 

1 r A i r A _ qA 2 

Vq = ^Lcfl(T = 0) = J dpp 2 (16\e\ + 16|e|) - J dpp 2 (2\e'\ + 2|e'|) + - G v p 2 , (18) 

with 


e = ±\/(p — p) 2 + A 2 , 


£ = ±^/(p + p) 2 + A 2 , 
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At this point, we observe that in the high-density region required for the existence of the CFL phase the leading 
parameter is the quark chemical potential p, whose contribution to the thermodynamic potential occurs only through 
diagrams with internal lines of fermions. This allows to integrate in the ghost fields and neglect the ghost contribution 
to the gluon polarization tensors. In a similar way, we can neglect all the gluon self-interactions in (fTTil) and dm). 
Moreover, the leading contribution in the infinity sum in (1121) comes from the term with the lowest order in the a s 
expansion. Therefore, the gluon effective action reduces to 


I dx^-^Gt - d v G*:f 



( 20 ) 


on which we used a covariant gauge F ,t = 0, with £ the gauge fixing parameter. 

The tensor in (TSUI) represents the fermion contribution to the one-loop gluon self-energy. In the hard-loop 
approximation j24| its leading part is fl5l ] 


n ^i?(Po = 0, p ->• 0) = [mf^o^o + rn 2 M S^iS^]S AB , (21) 

with the Debye and Meissner masses defined in Eq. m with the replacement p —> p,. Here, we introduced the 
notations rho = mr>0(A — p) + a/3 rh g 9(p — p)9{p— A) and tom = TOm$(A — p) to take into account the limitation to 
soft external momenta of lF£ B (jp) in the hard-loop approximation. Note that there is a different mass value depending 
on the momentum interval. For the Debye mass there are two values, % corresponding to the interval 0 < p < A, 
which is the infrared region of color superconductivity, and y/3m g corresponding to A < p < p, which is in the 
infrared domain of dense matter in the normal phase, where the Debye mass takes the known value y/3rh g f2~fl. |25|. 
We should note that, while the Debye mass is different from zero once there is finite density, the Meissner mass is 
only confined to the color superconducting domain p < A. 

With the help of EH . we can now integrate the gluons in (fTUl) to find 


8 -r-Fd 4 p E . * / ~ ~ m 

“G = sR indet[iA (p, m M ,m D )\, 


2/3 4- {2nf 


( 22 ) 


with gluon inverse propagator 
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and compact notation 



(24) 


for the three-momentum integral and sum in the Matsubara frequencies = = 7 jZL, n = 0,d=l,d=2,... The factor 8 in 
the right-hand side of (l22l) accounts for the degeneracy in the gluon spectra. 

Taking the Feynman gauge (£ = 1) and performing the Matsubara sum we obtain 


n G (r) = 



y/p 2 +m 2 D + 3\Jp 2 + m 2 M + /T 1 ln(l - e ^ p2+ ™°) + 3^ ln(l - e *V pa+ *") , 


(25) 


Therefore, in the T —> 0 limit, the gluon part of the thermodynamic potential becomes 

fl g = f2 G (T = 0) = T J dpp 2 [\Jp 2 + m 2 D 9(A - p) + 3rh 2 g 9{p - p)0{jp - A) + 3 \Jp 2 + m 2 M 9(A - p)^ (26) 

The thermodynamic potential of a gauged-NJL model has been also investigated in a recent paper [26|. However, 
this work did not consider color superconductivity and the study was done in the presence of a magnetic field. 

Combining all the contributions, we obtain that the zero-temperature, mean-field thermodynamic potential of model 
0 in the hard-loop approximation is given by 


fl - fig T fig 


H71 


where we subtracted the vacuum constant fl va c = fl(p = 0, A = 0). 

The dynamical quantities A and p should be found from the equations 


(27) 




(28) 


We call attention to two well-known but worth noticing features associated with the vector interactions of the 
model. First, as the mean value (ip”/o4’) enters in the covariant derivative as a shift to the particle chemical potential, 
the effective chemical potential for the baryon charge is now p instead of p. Second, while the solution of the gap 
equation (first equation in (1321) 1 is a minimum of the thermodynamic potential, the solution of the second equation is 
a maximum 0, since it defines, as usual in statistics, the particle number density p. 

Finally, let us define the model parameters that will be used in the numerical calculations of next sections. Following 
a standard procedure, we define the energy cutoff A = 602.3 MeV and the quark-antiquark coupling G 5 A 2 = 1.835 
and K A 5 = 12.36 to fit f„, m n , ttik and mg' to their empirical values in the sharp cutoff regularization [ 27 J. Then, 
the diquark coupling Gd , that produces a gap A ~ 10 MeV at p = 500 MeV, is found to be Gd = 1-2G 5 . A similar 
ratio Gd/Gs was already considered in [28f to investigate the M-R relationship in hybrid compact stars with color 
superconducting cores. Changing A in a few percentage, while simultaneously modifying Gn to produce the same 
value of A, does not affect our qualitative results. As for the values of the vector coupling, it is known that if the 
vector channel is originated from a Fierz transformation of a local color current-current interaction, the resulting 
coupling strength is Gy = 0.5 Gs- If instead, one starts from the molecular instanton liquid model or the PNJL 
model, the Fierz transformations give rise to values of Gy much more smaller [ 29 }. Based on these considerations, 
Gy is usually taken as a free parameter in the range Gy = (0 — 0.5)Gs. 


III. GLUONS EFFECTS ON THE EOS OF THE CFL PHASE OF STRANGE STARS 

In this section we shall explore the effects of the gluon rest-energy on the EOS of the CFL matter and the implications 
for strange stars. All our discussion will be based on the zero-temperature limit since compact stars typically have 
p T and our ultimate goal will be to investigate the star structure, which is insensitive to small temperature effects. 
The pressure P and energy density e can be found from the thermodynamic potential (l27l) as 

P = —(flq + fig — flvac) + ( B — Bo), e = flq + fig — fl va c + i^p ~ (B — Bq) (29) 

Notice that the chemical potential that multiplies the particle number density p in the energy density is p instead of 
p. This result can be derived following the same calculations of Ref. (3(j to find the quantum-statistical average of 
the energy-momentum tensor component Too, using for the present case the mean-field Lagrangian ©. 
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In (l29l) . we added, as usual, the bag constant B ( B 0 is introduced to ensure that e = P = 0 in vacuum). In the MIT 
bag model, B was introduced as a phenomenological input parameter to account for the free energy cost of quark 
matter relative to the confined vacuum (3l). Nevertheless, in the NJL model, the bag constant can be calculated in 
the mean-field approximation as a dynamical quantity related to the spontaneous breaking of chiral symmetry [32| . 
Following this dynamical approach, we have that the bag constant is given in our model by 


B = 





2 Gsiipi^i) 


+ 4 # (V’uV’uXV’dV’dX^aV’s) 


(30) 


where to,; is the dynamical quark masses, the quark condensates for the different flavors are given by 

n- , . 3 f A 2 , rn ? ; 

(V’iV’i) =-o / P dp- 


71-2 JpFi V m i + P 2 ’ 


(31) 


with ppi = ( 7 r 2 ^) 1 / 3 being the Fermi momenta depending on the densities of each flavor pi = (ipjipi); and the 
dynamical masses are found from the gap equations 


to, = 4 G s \ f p 2 dp + 2 

77 JpFi V m i+P 71 


where i ^ j ^ k take values in the flavor set u,d,s. 
The vacuum bag constant in (|29[) is given by 


p 2 dp 


rrii 


p 2 dp- 


m k 


yjm 2 + P 2 PFi y/rrii+p 


(32) 


B n = B 


Pu=Pd=Pa— 0 


(33) 


For the densities we are considering in the CFL phase (p > po), the phase transition to the restored chiral phase 
(with to, = 0) already took place. Hence, in our case B = 0 and the only remaining bag parameter in (1291) is Bq, 
given by 


Bq — —j 

T 


J p 2 dp ^ \Jm 2 +p 2 - \fp^ + 


2 Gsm 


y/ TO 2 + p 2 

with the same vacuum dynamical mass for the three quarks found from 


-4*1-2 


r d P p 2 - r ^= 

1 0 \Jm 2 +p 2 
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p 2 dp 


\/m 2 +p 2 
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-1 2 


p 2 dp 
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sj m 2 + p 2 


(34) 


(35) 


For the parameter set under consideration one obtains B 0 = 57.3 MeV/fm 3 [32 ]. 

We call attention that due to the lack of confinement in the NJL model, the contribution of the confined vacuum 
to the pressure and energy density cannot be consistently derived in the framework of this model, what always result 
in certain degree of ambiguity in the proposed value for the bag constant. 


A. Absolute stability 

The first thing we need to elucidate is what is the range of parameters that ensures absolute stability of the quark 
matter phase here considered. Absolute stability of the strange matter exists as long as the matter energy per baryon 
number at zero pressure 


E/A\p-q = 


(fig — flg ae + Hp) |p —Q + /IQ P 

P /3 


-fig + flg ae + (IQ P 
P /3 


(36) 


remains below the corresponding value for the iron nucleus (roughly 930 MeV). Here po is the value of the chemical 
potential at zero pressure. When the parameters are such that E/A < 930 MeV 1 the deconfinement of nuclear matter 
into this phase of strange matter is energetically favored. In (1361) we took into account the relation between the baryon 
and quark number densities ps = \{pu + Pd + p s ) = \p- 
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FIG. 1: (Color online) Energy per baryon number of quark matter at zero pressure as a function of the vector coupling. The 
horizontal long dashed line indicates the limiting value of 930 MeV (energy per baryon number of the iron nucleus). The solid 
line is for CFL matter without gluons, and the dashed line, for CFL matter with gluons. 


Fig. [I] shows the behavior of E/A\ P=0 with the vector coupling Gy for the CFL matter with and without gluons. 
The variation of the energy per baryon with the vector interaction at zero pressure is very smooth. This is expected 
because of the pressure normalization at low density, so the stability of the CFL matter is little affected by Gy. The 
effect of gluons on the energy per baryon number at zero pressure will depend on how much it lowers the pressure at 
a given density. If the effect were big enough to shift the zero pressure point to high values of quark density (high fi), 
the last term in Eq. (HhH) would push the energy up and destabilize the system. Nevertheless, the presence of gluons 
(dashed line) shifts the zero pressure point to a higher value of the density but not high enough as to take the system 
out of the absolute stability region. This shift being small, the zero pressure point is now closer to the minimum that 
usually exists in the energy vs. density curve, leading to a slightly lower energy per baryon number than without the 
gluons. 


B. EOS of strange quark matter with gluon contribution 

Fig. [2] displays the EOS for a range of Gy with and without the gluons contribution to the thermodynamic 
potential. The inclusion of gluons degrees of freedom soften the EOS, as can be seen by comparing the two panels. 
Only relatively high values of Gy can make the EOS of the CFL matter with gluons stiffer than the EOS of the 
regular CFL matter with no gluons and no vector coupling, at least in the energy density region relevant for compact 
stars. Of course, if we compare the curve for the CFL matter at certain value of Gy with the corresponding one 
including gluons and the same coupling value, the former will be stiffer than the second one, independent of the Gy 
value. 

We can also note how the shift in pressure due to a change in Gy for CFL matter with gluons for the same value of 
energy density increases with the increase of the energy density. While at energy densities smaller than 500 MeV/fm 3 
there is a negligible change even for Gy/Gs = 0.5, at e = 1500 MeV/frn 3 , the jump in pressure is really noticeable. 
This is an expected result, since the effect of the vector interaction is proportional to the baryon density and should 
be more significant in the high density region. 


C. Speed of sound in strange quark matter 

Another parameter that characterizes the dense matter is the speed of sound Vg = dP/de, which measures how 
fast a disturbance in the pressure travels in the dense medium. It can be used to constraint the EOS, since it has 
to obey some natural conditions like causality vs < c, and thermodynamic stability Vg > 0. A nonrelativistic gas 








FIG. 2: (Color online) Equation of state for CFL matter with gluon contribution (right panel) and without it (left panel), for 
different values of the vector coupling Gy. 




FIG. 3: (Color online) Sound speed as a function of the baryon density for CFL matter without the gluon term (left panel) 
and with it (right panel), for different values of the vector coupling strengths. The starting point to the left of each curve 
corresponds to the zero pressure point. The horizontal long dashed line determines the causality limit (vs = c). 


of nuclear matter has a very small speed of sound. The opposite extreme, an ultrarelativistic gas of free particles, 
has Vg/c 2 = and adding perturbative interactions or masses just leads to v’g/c 2 < \ am. Based on these facts, 
it has been speculated that for the range of intermediate densities meaningful for neutron stars, the speed of sound 
should lie somewhere between the values of these two extremes. These facts have led to conjecture the existence of a 
fundamental bound v 2 g/c 2 < | for the speed of sound of strongly interacting quark matter systems. However, as has 
been recently argued [34j], this bound seems hard to be reconciled with the observation of neutron stars with 2 Mq, 
at least for hybrid stars with all the reasonable low density EOS. 

In Fig. [3] the sound speed is plotted as a function of the baryon density for CFL matter without gluons (left panel) 
and with gluons (right panel). The sound speed is very sensitive to the presence of the vector interactions, as can be 
seen from the slopes of the curves in two panels of Fig. [3] Since the gluons have Debye and Meissner masses in the 
CFL phase, their presence tends to decrease the speed of sound for each given density and coupling Gy. As expected, 
a vector interaction tends to increase the speed since it hardens the EOS. Notice that despite the softening of the 
EOS by the gluons, the condition u|/c 2 < | is violated for all the vector interactions compatible with the absolute 
stability. 
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IV. GLUON EFFECTS ON THE M-R RELATIONSHIP OF STRANGE MATTER IN THE CFL PHASE 

An immediate application of the EOS just derived is to use it to construct stellar models of compact stars. In this 
regard, there are two distinct possibilities: the formation of strange stars with the strange quark matter in the CFL 
phase described by the EOS (OHl) . where the parameters are chosen to ensure absolute stability; or the formation of 
hybrid stars, on which the quark phase can be metastable, so the star cannot be all made of quark matter but it can 
exist in a core that is surrounded by nuclear matter. Given the results of Fig. |T]we are going to investigate the first 
possibility. With this aim, we should obtain the M-R sequence for strange stars made of matter with EOS given by 
(EH1) and explore how the gluons and vector interactions of the model can affect the sequence. Of particular relevance 
is to determine whether this sequence can reach the observed mass values of the compact objects, PSR J1614-2230 
and PSR J0348+0432 with M = 1.97 ± 0.04M© Q and M = 2.01 ± O.O4M 0 0 respectively. 

The M-R relationship is obtained by integrating the relativistic equations for stellar structure, that is, the well- 
known Tolman-Oppenheimer-Volkoff (TOV) and mass continuity equations, which in natural units, c = G = 1 are 
given by 


dM 

~dR 

dP 

dR 


47 tR 2 6 

eM / P\ / 47 tR 3 P\/ 2 M\~ l 

~1pV + ~){ 1 + ~W~) V - If) 


(37) 

(38) 


with P and e taken from (1291) . The corresponding M-R sequences are shown in Fig. [4] for EOS of quark matter 
with and without the gluon contribution. Comparing them, it is evident that the gluons decrease the maximum mass 
for each sequence up to 20%, this effect being more prominent for lower values of the vector interaction. Sequences 
including gluons within our model cannot reach 2 Mq if Gy/Gs < 0.2. Here, we should point out that several results 
actually suggest that a low vector coupling between quarks may be favored at high-densities [29] . Even more, as 
discussed in [35(, the vector interaction makes the chiral phase transition weaker in the low T and high p, region, with 
the possibility that the transition becomes a crossover in the region when the interaction is strong enough jl]]j Hsf, 
what is not expected from Lattice calculations. Then, the absence of the vector interaction would be preferable in 
the high density region. The same conclusion was reached by studying the phase transition from the hadronic phase 
to the quark phase in Ref. [37]. There, it was found that the hadron-quark phase transition may take place only 
at a small Gy/Gs ratio. For a larger ratio, the repulsive vector interaction makes the NJL phase too stiff to allow 
the crossing from the hadronic phase. However, such a phase transition could take place if one assumes zero vector 
interactions in the quark matter side and further considering effects that soften the EOS like diquark condensation. 

Although at present there is no definitive way to determine the value of the Gy coupling in the high-dense region, all 
these results provide an indication that Gy could be significantly small in the quark phase. If this is the indisputable 
case, then if Gy/Gs < 0-2, strange stars composed entirely of CFL quark matter with gluons, within the formalism 
presented here, should be disregarded, as they cannot explain the mass measurements of the mentioned compact stars 
PSR J1614-2230 and PSR J0348+0432. 


V. CONCLUSIONS 

In this paper we investigate the EOS of strange stars described by a three-flavor gauged-NJL model with vector 
and diquark channels. In the region of large densities and within the mean-field approximation, the ground state of 
the system is the same as in the CFL phase, but with the replacement of the baryon chemical potential by an effective 
chemical potential Jj. = p — Gyp. In addition, the EOS of the system contains the contribution of the rest energy of 
the gluons, which acquires Debye and Meissner masses in the color superconducting background through polarization 
operators with lines of fermions in the hard-loop approximation. The pressure and energy density defining the EOS 
were then used to calculate the M-R relationship of strange stars in such a parameter range. 

We have found that gluons soften the EOS of the CFL matter. This renders maximum stellar masses lower than 
the observed ~ 2 Mq if Gy/Gs < 0-2. This effect could exclude pure color superconducting strange stars as possible 
candidates of compact stars unless vector interactions were strong enough at high densities. New interactions, as the 
diquark-diquark repulsion studied in Ref. 38] or higher order quark interactions 39j, could bring enough outward 
pressure to make the favored composition of neutron stars to be pure superconducting quark matter with low vector 
interaction. Other quark-matter phases as 2SC with strong diquark coupling and including gluon contributions would 
have to be investigated as possible candidates for hybrid stars. 

We underline that while the energy eigenvalue e = 0 is always allowed in the partition function of non-relativistic 
bosons in a big box with periodic boundary conditions, so at zero temperature these particles always occupy the zero 
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FIG. 4: (Color online) Mass-radius relationship of strange stars made of CFL matter with (right panel) and without (left panel) 
the gluon contribution. The softening of the EOS by the gluons significantly decreases the maximum mass of the sequence, 
which for vector interaction with Gv/Gs < 0.2 remains below 2Mq. 


energy ground state with zero momentum, and hence zero pressure, things are very different for relativistic bosons. 
In the relativistic case, if the bosons have nonzero rest energy, they produce a positive contribution to the energy 
density and a negative term in the pressure. Since gluons in the color superconductor acquire rest energy due to their 
density-dependent Debye and Meissner masses m, these energies affect the EOS of the CFL matter at T = 0, thereby 
giving rise to all the interesting effects discussed in this paper. 

An important point to be indicated is that in the framework of the NJL model, the only vacuum contribution that 
can be calculated from first principles is that corresponding to the chiral condensation, as presented in Sec. Ill, while 
the additional contribution that will be related to confinement cannot be derived in a self consistent way. This is a 
limitation of the phenomenological NJL-type models used to describe the color superconducting stellar medium. 

To finish, we want to stress once again that due to the limitations of the used model, which is expressed in terms 
of a set of undetermined parameters (A, B, Gy, etc.), any statement regarding the existence or not of strange stars 
has to be taken more as an indication than as a definitive conclusion. 
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